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Abstract 


Using  an  impulse  or  Dirac  delta  approach,  the  stability  and  growth  of  yaw 
for  a  spin-stabilized  projectile  as  it  transits  a  yaw  card  range  is 
investigated.  The  card-induced  changes  in  the  complex  yaw  arm 
amplitudes  and  phases  are  expressed  as  difference  equations.  For  a  yaw 
card  range  with  uniform  spacings,  the  solutions  to  the  general  difference 
equations  yield  the  magnitude  and  phase  values  for  the  complex  yawing 
arms  as  a  function  of  the  card  index  number.  A  stepwise  encounter  v^th  a 
yaw  card  reduces  the  epicyclical  phase  values  across  the  yaw  card.  The 
parameters  for  the  solution  include  encounter  phase  value  at  the  entrance 
of  the  range,  the  characteristics  of  the  yaw  card,  and  the  distance  between 
the  cards.  Critical  curves  separating  stable  and  unstable  flight  are 
presented  as  a  function  of  yaw  card  spacing  and  a  stability  parameter  that 
depends  upon  both  the  flight  characteristics  of  the  projectile  and  the  card 
material  and  thickness.  Quasi-universal  curves  of  a  transformed  epicyclical 
phase  encounter  value  versus  a  normalized  epicyclical  phase  reduction 
value  are  graphically  presented. 
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YAW  CARD  INFLUENCE  ON  STABILITY  AND  YAW  GROWTH 
FOR  SPIN-STABILIZED  PROJECTILES 


1.  INTRODUCTION 

Yaw  card  arrays  set  up  at  equal  intervals  on  a  firing  range  have  been  used  to  obtain 
flight  data  when  spark  ranges  are  not  available  or  when  the  spark  range  operations  could 
be  compromised.  The  interactions  of  the  projectile  with  the  yaw  cards  cause  the  reduced 
flight  coefficients  to  differ  from  the  spark  range  values.  McCoy  (1992)  used  a  Fourier  analysis 
approach  and  assumed  that  the  fluctuating  details  generated  by  the  yaw  card  interactions 
would  average  zero,  thereby  obtaining  an  average  phase  shift  for  each  projectile-card  interac¬ 
tion.  These  projectile  flight  data,  when  reduced  and  corrected  for  the  yaw  card  interactions, 
predicted  the  flight  coefficient  values  obtained  in  a  spark  photography  range.  Cooper  and 
Fansler  (1997),  using  a  Dirac  delta  approach,  considered  the  individual  interactions  of  the 
projectile  with  each  yaw  card  and  calculated  the  fast  and  slow  arm  magnitudes  at  each 
card  along  with  the  phase  changes.  With  the  new  approach,  they  were  able  to  examine  the 
variations  from  the  average  value  at  each  card.  The  calculations  for  a  given  card  depended 
upon  the  calculation  obtained  for  the  previous  card.  In  turn,  the  previous  card  calculation 
depended  upon  the  adjacent  card  before  it.  Thus,  calculations  needed  to  be  performed  for  all 
the  cards,  beginning  at  the  first  card.  For  small  perturbations,  they  obtained  the  correction 
needed  to  reduce  the  yaw  card  range  flight  data  to  the  equivalent  spark  range  data,  which 
was  previously  obtained  using  McCoy’s  (1992)  approach.  Thus,  the  correction  is  only  valid 
for  small  perturbations,  or  equivalently,  sufficiently  thin  yaw  cards. 

Because  the  material  comprising  the  yaw  cards  is  dense,  it  is  reasonable  to  expect  that 
decreasing  the  spacing  of  the  yaw  cards  sufficiently  or  using  thick  yaw  cards  might  result 
in  destabilizing  the  flight  of  projectiles  through  the  range.  In  particular,  marginally  stable 
projectiles  would  require  careful  design  of  the  range  to  avoid  inducing  instability.  McCoy 
(1992)  asserted  that  the  usual  stability  criterion,  extended  with  his  correction  to  the  moment 
for  the  yaw  card  range,  was  still  valid.  Cooper  and  Fansler  (1997)  validated  the  assertion 
by  calculating  the  unstable  growth  that  met  McCoy’s  criterion  for  instability.  They  also 
pointed  out  that  the  obtained  expression  was  actually  an  approximation  valid  only  for  small 
values  of  the  stability  parameter,  which  is  equivalent  to  assuming  that  the  yaw  cards  are 
sufficiently  thin.  Cooper  and  Fansler  also  calculated  the  card-by-card  growth  of  a  projectile’s 
yaw  in  a  range  that  has  its  card  spacing  corresponding  to  nearly  the  optimum  condition  for 
maximum  amplitude  increase  for  the  modal  arm  amplitudes.  They  could  also  examine  the 
detailed  yawing  for  a  projectile  that  entered  the  range  at  different  initial  epicyclic  phase 


1 


values,  both  for  stable  and  unstable  flight. 

The  current  approach  is  a  continuation  of  the  work  by  Cooper  and  Pansier  (1997).  As 
before,  interactions  with  the  yaw  cards  are  treated  as  discrete  impulses,  but  now  a  yaw  card 
range  with  equal  intervals  between  yaw  cards  is  assumed.  This  approach  yields  second  order 
difference  equations  with  known  solutions.  With  solutions  for  the  second  order  difference 
equations  known,  it  is  no  longer  necessary  to  know  the  behavior  of  the  projectile  at  all  the 
preceding  yaw  cards.  In  this  report,  the  boundaries  between  stable  and  unstable  flight  are 
determined  along  with  the  growth  rate  for  yaw  modal  arm  amplitude  and  its  dependence 
on  various  parameters.  These  expressions  are  valid  for  both  small  and  large  values  of  the 
stability  parameter.  Quasi-universal  functions  are  obtained,  giving  the  relationships  between 
the  phase  encounter  value  at  each  card,  the  epicyclic  phase  jump  value,  and  the  change  in 
the  modal  arm  magnitudes  caused  by  the  card  encounter. 


2.  THE  DIFFERENTIAL  EQUATION  OF  MOTION  FOR  PROJECTILES 
FIRED  THROUGH  YAW  CARDS 


The  differential  equation  of  pitching  and  yawing  motion  for  a  spinning,  symmetric  pro¬ 
jectile,  acted  upon  by  a  linear  pitching  moment  and  neglecting  damping  processes  is  given 
in  Murphy’s  notation  (1963): 

e-iPi'-Mi=o,  (1) 

in  which 


C  = 

P  = 

M  = 
h  = 

ly  = 

d  = 
P  = 
P  = 
A  = 

Cmc  = 


sin  ^  +  i  sin  a,  the  complex  yaw, 

TyV 

PJ^r 

21  y 

axial  moment  of  inertia, 
transverse  moment  of  inertia, 
projectile  reference  diameter, 
axial  spin 
air  density, 

7rc?^/4,  reference  area, 

aerodynamic  pitching  moment  coefficient. 
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and  the  independent  variable  is  the  distance  in  calibers.  For  this  report,  it  is  assumed  that 
the  amplitude  of  the  complex  yaw  is  small  enough  that  sm  ct  ~  a  and  sin  /3  ^  /3.  Only  the 
first  order  aerodynamic  coefficients  are  then  used. 


Equation  (1)  is  the  equation  of  motion  in  air  but  the  yaw  card  range  has  cards  consisting  of 
dense  material  at  various  spacings.  Following  Cooper  and  Fansler  (1997),  a  nondimensional 
card-overturning  moment  can  be  defined  for  the  card  material,  similarly  to  M: 


(2) 


The  corresponding  nondimensional  impulse  for  the  card-overturning  moment  is  defined  as 


_  PcA(Ptc  ^ 

"  ~  2Iy 

and  the  independent  variable  is  the  distance,  s,  in  calibers. 


(3) 


For  the  yaw  card  range.  Equation  (1)  must  be  modified  to  account  for  the  yaw  card 
interactions.  As  discussed  in  the  preceding  section,  each  yaw  card  transmits  an  overturning 
moment  to  the  projectile,  which  can  be  modeled  as  an  impulse  located  at  the  position 
Sj,  which  is  the  arclength  traveled  by  the  projectile,  given  in  calibers.  The  Dirac  delta 
function  (Dennery  &:  Kryzwicki  1967)  is  used  with  the  pitching  moment  impulse  value,  and 
the  interactions  of  these  yaw  cards  with  a  projectile  are  represented  by  replacing  M  with 
M  -\-Xc  —  -Sj).  Equation  (1)  is  replaced  with 

r  -  iPi'  -  [M  J,  x:  S{s  -  s,)]e  =  0.  (4) 

j=i 


3.  DIFFERENCE  EQUATIONS  FOR  THE  DIRAC  DELTA  APPROACH 


The  solution  of  Equation  (4)  in  the  interval  between  the  jth  card  and  the  jth  -f  1  card 
is,  as  given  by  an  extension  of  Murphy’s  notation  (1963): 

i  =  -I-  j  =  0, 1,  •  •  • , 


,  n 


1, 


(5) 


in  which  (f>ij  and  <^2,j  are  the  card-modified  initial  phase  angles  of  the  fast  modal  arm  and  the 
slow  modal  arm,  respectively,  after  interaction  with  the  jth  yaw  card  but  before  interaction 
with  the  jth  -f-  1  yaw  card.  The  phase  angular  velocities  for  the  fast  and  slow  modal  arms, 
(f>\  and  <j)2,  are 

-  4M, 

=  if  -  -  iM. 
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(6) 

(7) 


It  is  assumed  that  the  projectile  is  spin  stabilized  and  is  gyroscopically  stable  while  traveling 
through  the  air.  This  stable  condition  requires  that  y/P'^  —  4M  is  real  and,  in  turn,  implies 
that  the  values  of  and  are  real.  To  achieve  stable  flight,  even  in  cold  and  dense 
air,  the  axial  angular  velocity,  p,  must  be  large  enough  so  that  a  comfortable  margin  of 
stability  is  achieved.  Good  margins  of  stability  are  achieved  for  most  projectiles  with  twist 
designs  varying  from  one  revolution  in  seven  bore  diameters  to  one  revolution  in  20  calibers. 
Likewise,  Kij  and  K2,j  are  the  fast  modal  arm  and  slow  modal  arm  amplitudes,  respectively, 
after  interaction  with  the  jth  yaw  card  but  before  interaction  with  the  jih  +  1  yaw  card. 
The  value  of  j  can  be  zero,  which  indicates  that  the  projectile  has  not  yet  interacted  with  a 
yaw  card.  Equation  (5)  describes  the  motion  of  an  epicycloid  in  which  both  modal  arms  are 
moving.  Equation  (5)  can  be  rearranged  to  illustrate  the  salient  aspects  of  the  motion  as  it 
relates  to  yaw  cards: 


J 


(8) 


in  which  we  use  the  relationship,  <f>[  —  <f>2  =  \/P^  —  4M.  Equation(8)  shows  that  if  the 
observer  moves  with  the  slow  arm,  a  circle  is  completed  in  the  distance  interval.  As  = 
(27r) —  4Af).  The  superimposed  circle  constitutes  an  epicycle.  For  a  statically  unstable 
projectile,  both  arms  move  in  the  same  direction.  The  phase  angle  value  of  the  first  term  in 
brackets  when  the  projectile  encounters  a  yaw  card  will  determine  the  impulsive  forces  on  the 
projectile.  Consider  when  ifij  =  K2,j.  When  the  epicyclic  phase  angle,  <l>i,j-(f>2,j+{<f>'i-<f>2)^^ 
equals  zero,  the  modal  arms  are  aligned  in  the  same  direction  and  d\(\/ds  =  0.  When  the 
phase  angle  equals  7r/2,  dj^j/ds  is  decreasing  most  rapidly.  When  the  phase  angle  equals  tt, 
the  total  yaw  and  d|||/ds  are  zero.  When  the  epicyclic  phase  angle  equals  3x/2,  d|^|/ds  is 
increasing  most  rapidly. 


For  convenience,  the  complex  modal  arms  are  defined  as 

Ku  =  ATwe'*'-, 

With  the  above  definitions.  Equation  (5)  becomes 


(9) 

(10) 

(11) 


Equation  (4)  can  be  integrated  across  the  jth  card  to  obtain  the  difference  equations: 

(12) 

=  Jd(si).  (13) 

The  negative  sign  denotes  the  side  of  the  yaw  card  that  the  projectile  approaches  and  the 
plus  sign  denotes  the  side  of  the  card  that  the  projectile  leaves.  Equations  (12)  and  (13) 
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show  that  the  value  of  the  epicyclic  yaw  is  not  changed  by  the  projectile-card  interaction, 
but  the  derivative  of  yaw  is  changed. 


3.1  Modal  Parameters  Before  and  After  Interaction 


Equation  (11)  and  its  derivative  with  respect  to  s  can  be  substituted  into  the  difference 
equations  and  solved  for  fCkj+i  in  terms  of  /Cfcj,  in  which  k  =  1,2  and  j  can  also  have  the 
value  zero: 


II 

c». 

1 

1 

ilc 

f>[-4 

(14) 

^2,j+l 

-  ^2,i  +  - 

ilc 

//  1! 

h  -  n 

(15) 

Cooper  &  Pansier  (1997)  defined 


Ic 


^  “  <I>1-4>'2 

Tc 

y/P2-iM' 

(16) 

-  ^'2)^3 

=  y/P'^  —  4:Msj, 

(17) 

so  that, 


(18) 

^2,3+1 

=  ^2,j  +  ifi}C2,j  +  5 

(19) 

in  which  fj,  is  called  the  impact  stability  parameter.  Equations  (18)  and  (19)  were  used  to 
obtain  the  point-by-point  results  and  the  approximation  that  delineated  the  line  between 
stable  and  unstable  flight  behavior  (Cooper  &  Pansier  1997). 


The  square  of  the  moduli  of  and  fC2,j+i  can  be  immediately  obtained  from  Equa¬ 

tions  (18)  and  (19). 

fe)  = 

/*ni  +  (^)"  +  2^cos(«,«)l,  (20) 

(tf)’  ■  ■--ft-w.'i* 

/[l  +  (|^)"  +  2|l^cos(*,„)).  (21) 
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=  5j+i  -f 

(22) 

^2,ji 

(23) 

in  which 

in  which 

and  is  called  the  epicyclic  phase  encounter  value  at  the  ;th  +  1  card. 

The  explicit  expression  for  the  sine  of  the  modal  arm  phase  angle  change  at  each  yaw 
card,  which  was  obtained  earlier  (Cooper  &  Fansler  1997),  is 

+  K2,jCos{<^j+i] 

-^1,1+1 

n[K2,j  +  Kij  cos($j+i)] 


sin(<^ij+i  —  <f>i,j)  — 

sm{<f>2,j+i  —  <f>2,i)  = 


K' 


2,j+l 


(24) 

(25) 


For  =  K2,j,  Equations  (24)  and  (25)  show  that  the  card-modified  phase  angle, 
of  the  fast  arm  is  always  reduced  by  the  impact,  while  the  card-modified  phase  angle,  (f>2j, 
of  the  slow  modal  arm  is  always  advanced.  The  total  epicyclic  phase  change  caused  by  the 
impact  would  be  ^j+i  —  The  epicyclic  phase  jump  is  defined  as 

(Tj+i  =  —  (<^j+i  —  ^j),  (26) 

Figure  1  shows  the  epicyclic  phase  jump,  and  the  modal  amplitude  change  induced  by 
the  interaction  normalized  by  their  maximum  values  as  a  function  of  the  epicyclic  encounter 
value,  $j+i. 


Figure  1.  Normalized  Phase  Jump  and  Modal  Arm  Amplitude  Versus  the  Epicyclic  Phase 
Encounter  Value  (^  =  0.1). 


The  maximum  value  of  o'j+i  can  be  easily  obtained  from  Equation  (25)  with  the  assump¬ 
tion  that  crj+i/2  =  4>2,i+i-<i>2,i-  Differentiation  of  Equation  (25)  with  respect  to  the  encounter 
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phase  value  and  setting  the  expression  equal  to  zero  will,  with  trigonometric  identities,  yield 
the  maximum  value,  <7j+i  =  4tan“^,  which  is  used  to  normalize  in  Figure  1. 

Also  superimposed  upon  the  figure  is  the  normalized  curve  for  the  fractional  increase  or 
decrease  in  the  modal  arm  amplitude.  The  phase  encounters  resulting  in  decreasing  modal 
arms,  <  0,  occur  for  2tan“^/<  <  <  tt,  while  the  phase  encounters  that  yield 

only  increasing  modal  arm  amplitudes  occur  for  the  values,  tt  <  <  27r  +  2tan“^/u. 

These  boundary  points,  or  points  where  Aifn  =  0,  can  immediately  be  obtained  by  setting 
the  amplitude  ratios  in  Equations  (20)  and  (21)  to  one  and  solving  for  the  phase  encounter 
values.  The  minimum  and  maximum  values  of  AAr„  occur  when  =  7r/2  +  tan~^  n  and 
$  =  35r/2  +  tan”^  fi,  respectively.  These  extreme  values  can  be  obtained  by  differentiating 
Equations  (20)  and  (21)  with  respect  to  the  phase  encounter  value  and  setting  equal  to  zero 
and  solving.  The  negative  part  of  the  curve  showing  the  fractional  decrease  in  the  modal  arm 
amplitudes  is  normalized  by  the  value,  n  —  \/l  +  /z^  +  l,  whereas  the  fractional  increase  in  the 
modal  arm  amplitudes  is  normalized  hy  +  \/l  +  —  1.  The  interaction  takes  place  tan“^  // 

past  the  points  where  the  projectile  is  pitching  most  rapidly  down  and  up,  respectively. 
The  interaction  yields  (Tj+i  =  2  tan~^  with  the  interaction  effectively  taking  place  around 
minimum  and  maximum  values  of  d|^|/ds.  The  force  that  occurs  on  the  projectile  at  these 
values  is  along  the  lateral  motion  of  the  projectile.  This  is  an  example  of  the  general  result 
that  a  given  impulse  changes  the  kinetic  energy  maximally  when  the  absolute  value  of  the 
dot  product  of  the  momentum  with  the  impulse  is  a  maximum. 


An  invariant  relation  between  the  magnitudes  of  the  two  modal  arms  is  applicable  to 
both  stable  and  unstable  yawing  motion.  Some  manipulation  of  Equations  (20)  and  (21) 
yields 


ti'2  7^2  _  jr^2 

~  ■^2,i+l  -  -  ^2,j- 


(27) 


The  above  equation  then  implies  that 


t^2  _  t^2  t^2 
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(28) 


which  is  a  hyperbola  curve  in  two-dimensional  space  for  the  two  modal  arms,  which  expresses 
an  invariant  quantity  for  the  interaction  process. 


3.2  Solution  for  Uniform  Yaw  Card  Range  With  Difference  Equations 

Equations  (20),  (21),  (24),  and  (25)  can  be  used  to  obtain  a  point-by-point  calculation 
of  the  modal  behavior  for  cards  with  any  arbitrary  spacing,  but  yaw  card  ranges  are  usually 
comprised  of  a  series  of  equally  spaced  cards  to  facilitate  obtaining  the  projectile’s  flight 
coefficients.  If  the  intervals  between  the  cards  are  all  equal,  the  phase  interval  characterizing 
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the  range  is 


(29) 


7o  —  Sj+i  Sj, 

a  relationship  between  the  yaw  arms  for  three  succeeding  cards  can  be  obtained: 

+  =  0,  (30) 

tC2j+2-2ge^'>^^]C2j+i  +  e^'^IC2j  =  0,  (31) 

in  which 

7  =  7o  -  2m7r,  (32) 

g  =  /isin(7/2)  +  cos(7/2).  (33) 

in  which  m  =  0  or  a  positive  integer  and  0  <  7  <  27r.  When  m  >  1,  the  cards  are  spaced 

so  far  apart  that  more  than  one  cycle  of  motion  occurs.  The  symbol,  7,  will  be  called  the 

principal  phase  interval.  The  use  of  the  principal  phase  interval  simplifies  the  analysis  and 
equations.  It  is  assumed  that  the  values  of  <^1,0  and  (f>2fl  are  the  initial  phase  values  of  the 
fast  and  slow  modal  arms,  respectively,  found  at  a  phase  interval,  70,  in  front  of  the  first 
card. 


To  solve  the  difference  Equations  (30)  and  (31),  define 

ic,,i  =  (34) 

(35) 

Equations  (30)  and  (31)  then  become 

Qi,j+2  ~  2gQi,j+i  +  Qi,j  =  0,  (36) 

Q2,j+2  —  2gQ2,j+i  +  Q2,j  =  0,  (37) 


which  are  in  a  standard  form  (Hildebrand  1968).  In  the  above  form.  Equations  (36)  and 
(37)  are  equivalent  and  have  the  same  general  solution.  Another  standard  form,  which  may 
be  more  convenient  to  work  with,  is  provided  by  substitution  with  the  definitions; 

cos  6  =  g  for  g  <  1,  (38) 

cosh^  =  5'for^>l.  (39) 


When  ^  <  1,  the  difference  equations  yield  stable  solutions,  whereas  when  g  >1,  unstable 
solutions  are  obtained  (Hildebrand  1968).  For  S'  <  1,  the  general  solution  to  Equations  (36) 
and  (37)  is,  using  the  above  definitions, 

Qi,j  =  Ai,i  sin(j  -  1)^  +  Ai,2sin 
Q2,j  =  ^2,1  sin(i  -  1)^  +  ^2,2  sin  j0. 
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(40) 

(41) 


For  >  1,  substitute  iQ  for  9  to  obtain  the  hyperbolic  functions. 


Substitution  of  Equations  (40)  and  (41)  into  Equations  (36)  and  (37)  with  i  =  0  and 
j  =  1  gives  the  expression  for  the  modal  complex  yaw  arm  and  yields  for  ^  <  1,  in  terms  of 
the  initial  conditions: 

— {[^i,o(l  -  */^)  -  sin  sin(y  -  1)0},  (42) 

Sint/ 

e«7(i-i)/2 

^2,i  =  — r-^{[Xi2,o(l  +  iy)  +  sin;0  -  /C2,oe‘'^'^^  sin(i  -  1)0},  (43) 

sm  u 

The  quantities  within  the  square  brackets  in  the  above  equations  are  the  values  of  the  modal 
arms  after  passing  through  the  first  card.  Again,  for  y  >  1,  the  above  trigonometric  functions 
are  replaced  by  the  respective  hyperbolic  functions: 

-t7(j-l)/2 

=  RiTib0  “  e>/C2,oe"’^]  sinh  j0  -  sinh(j  -  1)0},  (44) 

e«7(i-i)/2 

^2,i  =  sinh0  ■*"  sinh  j0  -  sinh(;  -  1)0},  (45) 


The  ratios  of  the  jth  yaw  arm  moduli  to  their  original  values  can  be  obtained  either  from 
first  principles  or  with  algebraic  manipulation  of  Equations  (42)  and  (43),  with  substitution 
for  trigonometric  identities.  The  resulting  equations  are,  foTg<  1, 

^  ^5^[Cisinj0  +  Asin(i  -  1)0],  (46) 

\Ai,o/  Sin  0 

(^)  ^  J!^[C2sin;0  +  I>2sin(i  -  1)0],  (47) 

in  which,  the  constants,  Ci,  C2,  T>i,  and  7)2  are 


Cl 

C2 

A 

V2 


-2/1^^^  sin  $1  + l  +  ( 

jti.o 

-2)U^^^sin$i  + l  +  ( 
■ft.  2,0 

2/^^sin($i  -7/2), 

^1,0 

2/^^sin(6i  -7/2), 

^2,0 


^)2  +  2^cosi. 

^1,0  ^1,0 

|!£)^  +  2^cos4. 

^2,0  A  2,0 


(48) 

(49) 

(50) 

(51) 


in  which 

1*1  =  7  +  -  <f>2,o-  (52) 

Since  $1  corresponds  to  the  initial  conditions  at  the  first  card,  the  initial  starting  condition 

A 

can  be  relaxed  to  any  position  in  front  of  the  first  card  and  $1  can  be  replaced  by  $1,  as  in 
Equation  (22),  in  which 


=  *?1  +  <^0- 


(53) 


From  Equations  (20)  and  (21),  the  values  C\  and  C2  are  recognized  as  the  square  of  the 
moduli  minus  1. 


For  unstable  yawing  motion,  or  equivalently,  g  >  substitute  the  hyperbolic  functions  for 
the  trigonometric  functions  wherever  the  function  argument  is  0.  With  further  manipulation, 
the  expression  can  be  written  as  an  increasing  term  plus  a  decreasing  term  minus  a  constant 
value  term: 


K 


K 


2,j 


-  1 


(Ci  +  Pie-^)  +  (Ci  +  Dje^)  -  2  (Ci  +  A  cosh^) 

4  sinh^  9 

[C2  +  P2e-^)  +  [C2  +  -  2  (C2  +  V2  cosh  9) 

4  sinh^  9 


(54) 


(55) 


For  =  1,  the  general  solution  for  the  squared  amplitude  of  the  modal  arms  can  be  obtained 
by  taking  the  limit  of  9  going  to  zero  in  Equations  (46)  and  (47): 

2 

I  -1  =  i[(Ci  +  Di)i-I>i],  (56) 


'K,,i 


K- 


2J 


,^2,0. 


-  1 


i[(^2  +  7?2)i  ■“ 


(57) 


Although  the  projectile  flight  is  unstable  for  =  1,  the  modal  arm  amplitudes  increase  only 
linearly  instead  of  exponentially,  as  occurs  for  >  1. 


For  stable  yawing  motion.  Cooper  and  Fansler  (1997)  observed  that  the  values  for  the 
change  from  the  initial  modal  lengths  lay  in  an  envelope  with  a  midpoint  average  value  that 
was,  in  general,  not  zero.  If  Equations  (46)  and  (47)  are  integrated  with  j  treated  as  a 
variable  and  then  divided  by  the  integrated  interval,  the  average  values  are  found  to  be 


'K 


1  = 


1 


K-. 


2,i 


.^2,0, 


1  = 


2  sin^  9 

1 

2  sin^  9 


[Cl  “I"  7^1  cos  0\ , 


[C2  “I"  7^2  cos  0\ 


(58) 

(59) 


Here  the  average  value  depends  both  on  the  initial  modal  phase  and  the  card  interval  spacing. 
As  7  is  varied  so  that  the  stability  boundary  is  approached,  then  sin  9  becomes  smaller  and 
the  absolute  values  of  the  average  quantities  can  become  larger. 


4.  STABILITY  BOUNDARIES  AND  GROWTH  IN  YAW 

In  the  last  section,  it  was  shown  that  the  boundary  for  stability  corresponded  to  the 
condition  that  g  =  1.  With  trigonometric  identities,  it  can  be  established  from  the  definition 
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for  g,  Equation  (33),  that  the  stability  boundaries  for  7o  encompassing  phase  difference 
intervals  are  given  by 

Tiow  =  2m7r,  (60) 

Thigh  =  2m7r  +  4tan“^ /u,  (61) 

in  which  m  =  0  or  a  positive  integer  and  Kijq  =  K^fi-  The  usual  approximate  launch 
condition  from  an  unworn  gun  is  A'i,o/i^2,o  =  !•  The  subscripts  “low”  and  “high”  designate 
the  minimum  and  maximum  values  of  the  phase  interval  for  a  given  region  designated  by 
the  integer,  m. 

The  maximum  growth  curves  in  the  //,  70-plane  are  obtained  by  differentiation  of  Equa¬ 
tion  (33)  with  respect  to  7  and  setting  the  resulting  expression  equal  to  zero: 

7max  =  2m7r  -t-  2  tan“^  pi.  (62) 

in  which  =  -^"2,0-  Figure  2  shows  the  first  two  of  the  regions  for  stable  and  unstable 
behavior.  The  growth  rate  is  a  maximum  when  To  =  2  tan“^  pi  for  the  lower  unstable  region 
and  To  =  27r  -|-  2  tan"^  pi  for  the  next  unstable  region.  This  maximum  growth  rate  curve 
corresponds  to  =  y/\  -|-  pi^. 


Figure  2.  Stable  and  Unstable  Regions  as  a  Function  of  pi. 

Cooper  and  Fansler  (1997)  have  already  shown  that  with  some  assumptions  including 
//  <C  1,  the  condition  for  stability  as  given  by  McCoy  (1992)  can  be  reproduced.  This 
condition  can  also  be  obtained  with  the  present  approach,  assuming  only  that  //■Cl.  The 
conditions  for  instability  in  the  present  approach  are  that  ^  >  1  or,  equivalently, 

//sin  ^ -h  cos  ^  >  1.  (63) 
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Expand  the  trigonometric  functions  in  series  for  small  7.  For  limiting  small  values  of  fi,  only 
the  first  terms  need  to  be  kept  to  obtain  the  conditions  for  instability: 

0.  (64) 

With  Equations  (29),  (64),  (16),  (3),  and  (2),  and  some  algebraic  manipulation,  the  approx¬ 
imate  gyroscopic  stability  factor  for  stable  yawing  motion  is 

Sg  =  >  1-  (65) 

4{M  +  McTe/dc) 

The  boundaries  obtained  for  the  approximation.  Equation  (64),  and  the  exact  solution 
are  shown  in  Figure  3.  Here,  it  is  shown  that  the  approximate  solution  agrees  well  for  fx  <  0.4 
and,  for  conventional  yaw  card  ranges  [jj  =  0.1  for  a  card  range  designed  for  a  7.62- mm  rifle, 
according  to  McCoy  [1992]),  the  approximate  solution  is  adequate.  For  larger  values  of  fx, 
the  straight  line  approximation  of  7  by  Equation  (64)  is  larger  than  the  exact  value. 


Figure  3.  Comparison  Between  Approximate  and  Exact  Solution  for  the  Stability  Boundary. 

Some  general  conclusions  can  be  obtained  by  examining  the  unstable  behavior  of  the 
modal  arms  for  a  large  number  of  yaw  cards.  To  investigate  the  asymptotic  behavior  of 
card  interactions,  treat  the  card  index  number,  j,  as  a  continuous  variable  and  differentiate 
Equation  (44)  with  respect  to  j  obtaining 

^  (66) 

dj  2 

Integrating  Equation  (66),  it  is  obtained  that  for  a  card  index  value,  k,  larger  than  j  the 


12 


value  for  the  first  modal  arm  is  approximately 

(67) 

Equation(67)  shows  that,  for  large  index  card  numbers,  the  magnitude  of  the  modal  arms 
will  always  increase  and  the  phase  change  at  each  encounter  is  simply  half  the  principal  value 
of  the  phase  difference  interval.  Likewise,  the  second  modal  arm  amplitude  will  increase  in 
the  same  way,  but  the  second  modal  arm  will  be  advanced  instead  of  retarded.  These  results 
show  that  the  behavior  of  the  modal  arms  at  large  values  of  j  is  independent  of  the  initial 
conditions.  Note  that  no  restrictions  have  been  placed  on  the  relative  phases  or  sizes  of  the 
modal  arms.  Even  though  the  initial  phase  encounter  value  may  be  such  that  the  amplitudes 
of  the  modal  yaw  arms  decline  as  the  first  part  of  the  yaw  card  range  is  traversed,  the  phase 
encounter  values  eventually  adjust  so  that  the  phase  shift  and  amplification  factor  depend 
only  on  the  stability  parameter,  /x,  and  the  principal  value  of  the  phase  difference  interval, 
7.  For  small  values  of  //,  6'^  «  7(/x  —  7/4)  and  2(sr  -  1). 

For  large  values  of  j,  A'lj+i  ~  7(^2, j+i ,  from  Equation  (28).  With  Kij+i/Kij  approaching 
a  constant  value,  according  to  Equation  (66),  the  epicyclic  phase  encounter  value,  $j+i,  also 
approaches  an  unchanged  or  stationary  value,  from  Equations  (20),  (21),  (24),  and  (25).  This 
means  that  the  epicyclic  phase  jump,  cr^+i,  asymptotically  approaches  the  phase  interval 
value,  7,  as  the  projectile  advances  through  the  range.  The  value  of  would  depend  on 
the  value  of  7  «  determined  by  the  functional  form  shown  in  Figure  1. 

The  same  functional  form  is  obtained  from  the  solution  to  the  difference  equation.  Setting 
the  second  term  in  parentheses  on  the  right-hand  side  in  Equations  (54)  and  (55)  equal  to 
zero  and  solving  for  the  relationship  between  $1  and  7  gives  the  condition  for  stationary 
phase  encounter  values,  =  $1,  that  yield  modal  arm  amplitudes  with  a  constant  factor 
for  the  increase.  In  the  same  manner,  setting  the  first  term  in  parentheses  on  the  right-hand 
side  in  Equations  (54)  and  (55)  equal  to  zero  and  solving  for  the  relationship  between 
and  7  gives  the  condition  for  stationary  phase  encounter  values  that  yield  decreasing  modal 
arm  amplitudes.  The  value  of  $1  that  satisfies  these  conditions  for  a  given  value  of  7  is 
designated  as  One  value  of  would  be  found  for  both  an  increasing  and  a  decreasing 
value  for  the  change  in  the  modal  arm  amplitude. 

The  solutions  for  j„  are  given  in  Figure  4,  which  shows  the  phase  interval  values  as  a 
function  of  for  /x  =  0.1.  Equation  (67)  shows  that  the  phase  encounter  values,  2tan“^  /x  < 
<  X,  cannot  exist  for  large  values  of  j  with  cards  uniformly  spaced.  The  curve  is  shown 
as  dashed  wherever  there  is  a  fractional  decrease.  Because  the  phase  encounter  values  cannot 
occur  for  large  values  of  j  when  the  amplitude  is  decreasing,  these  corresponding  phase 
encounter  values  are  more  accurately  designated  as  pseudo-stationary  phase  solutions  or  the 
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phase  encounter  solutions  for  the  pseudo-stationary  branch  for  the  phase  interval  between 
yaw  cards,  7„. 


Figure  4.  Stationary  Phase  Encounter  Values  Versus  the  Phase  Interval  with  the  Superim¬ 
posed  Curve  for  the  Fractional  Change  in  the  Modal  Arm  Amplitudes  (/x  =  0.1). 

Much  of  the  information  about  stationary  phase  solutions  could  also  be  obtained  from  the 
single  card  treatment  of  the  preceding  section.  If  one  assumes  that  7  =  aj+i,  Figure  4 
can  be  obtained,  but  with  no  immediate  knowledge  of  the  pseudo-stationary  solutions.  The 
value  of  7  is  normalized  by  dividing  by  its  maximum  value,  4tan”V)  which  occurs  at 
=  2tan“^  n-  Figure  4  shows  that  the  modal  arm  amplitude  remains  constant  at  both 
=  2tan~^/x  and  =  tt.  The  latter  value  corresponds  to  the  position  where  the 
projectile  is  aligned  with  the  flight  direction  (neglecting  the  yaw  of  repose),  and  minimal 
interaction  with  the  target  would  be  expected. 

The  expression  in  the  second  set  of  parentheses  in  Equation  (54),  when  set  equal  to 
zero  to  obtain  stationary  encounter  phase  solutions,  can  be  used  to  obtain  a  relationship 
between  a  shifted  value  of  the  encounter  phase  value  and  the  phase  interval  between  cards. 
The  expression  in  the  first  set  of  parentheses  could  also  be  used  because  the  modal  arms 
are  assumed  equal,  and  a  squaring  operation  will  be  used  that  enables  applicability  to  both 
decreasing  and  increasing  values  of  the  modal  arms.  When  the  resulting  equation  is  expressed 
in  terms  of  —  7/2  and  the  results  squared  to  obtain  a  quadratic  equation  that  is  solved, 
a  transcendental  relationship  is  obtained  between  the  transformed  value  of  the  encounter 
phase  and  phase  interval  between  cards: 

cos$.y  =  _  cos(7/2).  (68) 

/i 

With  the  substitution  of  tan  o  for  n  and  using  trigonometric  identities.  Equation  (68)  can 
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be  represented  as 


cos 


sin(7/2  —  a) 


sin  a 


(69) 


in  which  a  =  tan“^  fi.  Both  the  quantities,  7/2  —  a  and  a,  in  Equation  (69)  are  small  enough 
so  that  they  can  well  approximate  their  sine  functions.  Finally,  dividing  7  by  the  normalizing 
factor,  4tan“^  /i,  yields  the  quasi-universal  relationship  between  the  normalized  value  of  7 
and 

7„  «  (1 +  cos$i)/2,  (70) 


in  which  —  ^/2  which  applies  for  all  realistic  values  of  fj.  for  a  yaw  card  range.  By 

transforming  to  the  variable,  the  relationship  between  the  phase  interval  difference,  7, 
for  the  card  spacing  and  the  value  of  needed  to  maintain  a  constant  value  of  the  epicyclic 
phase  encounter  value,  $j+i,  at  the  arbitrary  jth  +  1  card  is  given  by  Equation  (70). 


Comparison  with  Equation  (68)  for  the  exact  value  and  the  approximate  expression. 
Equation  (70),  shows  that  the  error  in  the  approximation  is  less  than  0.3%  for  fi  =  0.3  and 
for  the  transformed  values  at  Ri  rmr/A,  where  the  most  error  occurs  with  odd  integer 
values  of  m.  By  expansion  into  series  form,  it  can  be  shown  that  the  error  condition  when 
using  the  above  approximation  is 


Error  <  //^sin(2$i)(sin$i)/24.  (71) 

The  approximation.  Equation  (70),  becomes  exact  where  $1  =  m?r/2,  in  which  m  equals  an 
integer  with  zero  included. 

Equation  (70)  has  far-reaching  consequences  for  isolated  encounters  of  a  projectile  with  a 
yaw  card.  Since  CTj+j  =  7  when  the  condition  of  a  constant  phase  encounter  value  is  achieved 
and  the  card  encounter  phase,  determines  the  value  of  phase  encounter  jump,  crj^i,  it 

is  deduced  that 

CTj+i  2(tan“^  /^)(1  +  cos  $j+i),  (72) 

in  which  —  <Tj.,.i/2  for  an  arbitrary  yaw  card.  The  transformed  phase  encounter 

value,  may  be  interpreted  as  the  center  point  of  the  encounter  phase  value  around 

which  the  interaction  takes  place.  Immediately  after  the  interaction,  the  phase  encounter 
value  is  reduced  to  —  o-j+i . 

Figure  5  shows  the  normalized  value  of  the  phase  jump,  (T„,  for  /x  =  0.05  and  fx  =  0.3, 
plotted  as  a  function  of  the  transformed  phase  encounter  value,  $j+i.  These  values  of  the 
stability  parameter  range  from  one  of  the  smaller  values  to  one  of  the  larger  values  of  fx 
that  might  be  encountered.  The  special  case  for  stationary  caxd  encounter  phase  through  a 
uniformly  spaced  range  is  also  described  by  the  results  of  Figure  5,  in  which  7^  =  cr„.  The 
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differences  in  the  values  of  cr„  plotted  for  the  two  different  values  of  the  stability  parameter, 
H,  cannot  be  seen. 


Figure  5.  Normalized  Phase  Jump  and  Modal  Amplitude  Change  Versus  Quasi-Universal 
Transformed  Card  Encounter  Phase  Value. 

Also  shown  is  the  fractional  change  in  the  modal  arm  amplitude,  which  is  plotted  for  the  sta¬ 
bility  parameters,  //  =  0.05  and  fi  =  0.3.  The  plotted  curves  for  ^Kn  show  more  differences 
than  the  curves  for  7„.  These  curves  can  be  directly  plotted  by  expressing  Equation  (20)  in 
terms  of  the  variable, 

^ 

Equation  (72)  not  only  represents  a  simple  physical  basis  for  the  encounter  phase  and  the 
epicyclic  phase  jump  but  can  also  be  used,  because  of  its  fundamental  functional  form,  to 
derive  a  simple  expression  for  the  phase  encounter  value,  ,  and  to  show  how  it  approaches 
the  stationary  phase  encounter  value,  which  depends  on  7.  For  nonstationary  phase 
encounter  values,  is  simply  the  sum  of  the  encounter  phase  at  the  jth  -|-  1  card  plus 
the  interval  pheise,  7,  traveled  minus  the  phase  jump,  (Tj+i,  from  encountering  the  jth  -|-  1 
card: 

$i-|-2  =  ^j+1  +  7  ~  ^j+i- 

The  behavior  of  the  transformed  phase  encounter  value  with  card  index  number  can  be 
obtained  by  first  substituting  in  the  right-hand  side  of  Equation  (74)  with  the  results  of 
Equations  (70)  and  (72)  to  obtain 

—  ^j+i  =  2(tan~^  /x)(cos  —  cos  $j+i).  (75) 

If  the  phase  encounter  value,  corresponds  to  a  value  in  the  psuedo-stationary  region, 

Equation  (75)  dictates  that  the  value  of  the  next  encounter  phase  moves  away  from  the  value 
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of  the  encounter  phase  on  the  pseudo-stationary  branch  corresponding  to  7  and  toward  the 
encounter  phase  corresponding  to  7  on  the  stationary  branch.  Thus,  is  the  solution  for 
the  stationary  phase  encounter  value,  which  is  restricted  to  tt  <  <  27r.  Now  for  small 

values  of  //,  the  approximation,  —  ^j+i  ~  ^j+2  —  and  treat  j  as  a  continuous 

variable  to  obtain  a  differential  equation  whose  solution  is 


$j4.i  =  2  tan  ^ 


$ 

Itan^l 

*■ 

tany- 

tan  1  tanh( j  |  sin  $1 

L 

tan  ^  n) 

$ 

1  tan  1  —  tan 

Li 

.  —  tanh(y  1  sin  $1 1  tan' 

■V) 

(76) 


in  which  $1  is  the  initial  value  for  the  transformed  encounter  phase  and  j  assumes  the  value 
zero  and  any  positive  integer.  Except  for  the  singular  point  at  7r/2,  Equation  (76)  shows 
that  the  asymptotic  value  for  the  transformed  phase  encounter  value  is  always  negative, 
which  is  equivalent  to  the  asymptotic  phase  value  being  between  tt  and  27r.  These  results 
axe  consistent  with  Equation  (67),  which  shows  that  the  phase  jump  approaches  a  constant 
value  and  that  >  1  for  large  values  of  j. 


The  epicyclic  phase  at  card  encounter  as  a  function  of  the  card  index  number  is  shown 
in  Figure  6.  The  epicyclic  phase  at  encounter  with  a  yaw  card  was  calculated  point  by 
point  from  the  first  card  onward  (Cooper  &  Pansier  1997)  and  also  calculated  with  the  ap¬ 
proximation  solution  given  by  Equation  (76).  Comparisons  between  exact  and  approximate 
calculations  showed  no  significant  differences.  The  foregoing  discussion  describes  the  phase 
encounter  values  when  the  modal  arms  are  equal.  An  expression  has  been  obtained  for 
unequal  arms  but  is  not  given  here  because  of  its  extreme  complexity,  which  limits  insight. 


Figure  6.  Phase  Encounter  Value  Versus  Card  Index  Number. 
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5.  FURTHER  DISCUSSION  AND  ILLUSTRATIONS 


The  equations  developed  are  used  to  examine  the  detailed  yawing  motion  of  the  projectile 
through  a  yaw  card  range.  The  equations  giving  the  yaw  arm  values  for  each  interval, 
Equations  (42)  and  (43),  are  substituted  into  Equation  (5)  and  plotted.  The  first  example 
shows  the  yawing  motion  for  a  7.62-mm  projectile  passing  through  a  dense  yaw  card  range 
discussed  previously  (McCoy  1992;  Cooper  &  Pansier  1997).  Figure  7  shows  the  yawing 
motion  comparison  of  a  projectile  passing  through  10  cards  with  one  set  of  cards  having 
/i  =  0.1  while  the  other  set  has  another  projectile  passing  through  cards  having  n  =  0.0. 
Having  cards  with  (x  =  0.0  is  an  artifice  that  enables  tracking  the  trajectory  even  though  no 
cards  exist  in  the  range. 


a/K2_o  k,,o/K2,o-i 


Figure  7.  Comparisons  of  Stable  Epicyclic  Motion  for  the  7.62-mm  Projectile. 

The  fast  arm  and  slow  arm  are  initially  equal  in  length,  which  is  approximately  the  initial 
conditions  for  many  projectiles.  The  curves  begin  at  the  origin  with  =  27r  and  <^2,0  =  tt, 
which  yields  for  the  initial  value  for  the  epicyclic  yaw,  ^0  =  tt.  The  projectiles  travel  a  phase 
distance.  Si,  equal  to  the  phase  interval,  7  =  1.65,  so  that  $1  =  1.65-|-7r.  The  phase  encounter 
value  at  the  first  card  plus  the  parameter  values  given  in  the  figure  completely  define  the 
motion.  After  each  impact,  the  curves  alternate  between  dashing  and  continuous  lines.  The 
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curves  pass  through  the  origin  on  each  epicycle,  because  Kij  =  K2J  from  Equation  (28) 
when  Kifi  =  K2fl- 

Figure  8  shows  the  fast  axm  modal  amplitude  increasing  by  the  maximum  amount  with 
each  card  encounter.  The  projectile  encounters  every  card  with  the  phase  encounter  value 
$1  =  (4m  +  3)7r/2  +  tan“^ //.  For  this  phase  encounter  value,  cTj+i  =  2tan'"^/z,  and  the 
principal  phase  interval  must  be  such  that  7  =  2  tan~^  [x  to  obtain  maximal  amplification. 
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Figure  8.  Maximum  Growth  of  Fast  Modal  Arm. 


Figure  9  shows  the  fast  arm  modal  amplitude  changing  with  maximal  decay  for  each 
card  encounter.  The  projectile  encounters  every  card  with  the  phase  value  =  (4m  + 
1)77/2  +  tan“^  p.  As  discussed  before,  because  <Tj+i  =  2tan“^p  for  each  encounter,  the 
spacing  for  the  principal  phase  interval  must  be  such  that  7  =  2  tan“^  p  to  obtain  maximal 
diminishment  of  the  fast  modal  arm  amplitude.  Thus,  the  same  range  can  be  used  to  obtain 
Figures  8  and  9.  The  modal  arm  decrease  is  shown  for  only  12  cards  in  Figure  9,  but  with 
the  computer  precision  used  for  the  calculation,  the  curve  will  start  ascending  if  a  hundred 
or  more  cards  are  traversed.  The  first  terms  in  Equations  (54)  and  (55)  become  large  relative 
to  the  second  terms  for  large  enough  values  of  j.  The  smaller  the  precision  in  the  computer’s 
approximation  for  the  starting  value  or  the  principal  phase  interval,  the  sooner  the  curve 
ceases  descending  and  starts  evolving  into  an  exponentially  increasing  curve.  This  behavior 
has  its  counterpart  in  solutions  of  linear  differential  equations  with  exponential  solutions 
where  the  positive  exponential  term  has  a  vanishing  small  initial  component,  which  grows 
large  for  large  values  of  the  independent  variable.  In  practice,  it  would  be  difficult  to  design 
an  experiment  that  would  yield  the  results  shown  in  Figure  9.  Substantial  variations  from 
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the  initial  value  would  result  in  a  growing  amplitude  for  the  modal  arms  after  traversing  a 
few  cards. 

A  card  range  designed  to  produce  instability  could  be  useful.  For  instance,  high  values  of 
yaw  could  be  induced  by  judicious  placement  of  the  yaw  cards.  The  encounter  phase  values 
for  the  yaw  cards  can  also  be  controlled  with  the  yaw  card  range  design.  Cameras  could  be 
set  up  near  the  range  exit  to  study  the  flow  around  the  projectile  at  high  angle  of  attack. 


Figure  10  shows  the  detailed  epicyclic  yawing  motion  when  the  interval  corresponds  to 
the  maximum  value  of  g,  which  gives  the  maximum  growth  in  yaw  amplitude  through  a 
card  range  when  the  first  encounter  corresponds  to  the  maximum  yaw  amplitude  increase. 
Here,  the  starting  condition  is  <^o  =  and  after  traveling  somewhat  more  than  one  epicycle, 
encounters  the  first  card  at  =  Sx  +  2  tan“^  g.  The  yaw  amplification  at  the  first  card  will 
be  very  small  because  this  phase  value  for  the  first  encounter  is  very  near  to  the  encounter 
phase  value  for  no  amplification,  =  Stt.  Equations  (24)  and  (25)  show  that  the  delay  in 
the  total  phase  value  caused  by  the  encounter  with  the  first  card  is  very  small  for  /i  =  0.1. 
The  interval  value  is  larger  than  2x  and  in  subsequent  card  encounters,  the  encounter  phase 
changes  toward  the  value  =  2^^  +  37r/2,  which  corresponds  to  the  maximum  yaw 

amplification  increase.  The  yawing  motion  passes  through  the  origin  because  Ki^  =  K2,q. 
If  the  epicyclic  arms  had  been  unequal  initially,  the  ratio  of  the  arm  lengths  would  have 
approached  unity  as  the  yaw  card  range  was  traversed.  Then  the  point  of  closest  approach 
to  the  origin  for  each  epicycle  would  have  been  successively  less  as  the  yaw  cards  were 
traversed.  For  =  27r,  the  succeeding  encounter  phase  values  will  again  relax  to  the  phase 
value  that  gives  the  maximum  amplification.  For  any  given  initial  phase  value,  with  the 
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exception  of  =  7r/2,  the  succeeding  phase  encounter  values  will  relax  toward  the  phase 
value  that  gives  the  maximum  amplitude  increase. 

a  /  K2  0 


p  /  K2  0 


Figure  10.  Unstable  Epicyclic  Motion  for  a  7.62-mm  Projectile. 

6.  SUMMARY  AND  CONCLUSIONS 

Stability  and  yaw  growth  for  spin-stabilized  projectiles  passing  through  an  array  of 
equally  spaced  yaw  cards  are  investigated.  The  encounters  of  the  projectiles  with  the  cards 
are  treated  as  impulses,  or  mathematically,  as  Dirac  delta  functions.  Solutions  obtained 
by  considering  only  one  card  yield  expressions  for  modal  arm  amplitude  changes  and  also 
phase  changes  that  depend  only  upon  the  yaw  card  encounter  phase  value  and  the  stability 
parameter. 

The  encounter  with  the  yaw  card  retards  the  fast  arm  phase  and  accelerates  the  slow 
arm  phase,  resulting  in  a  reduction  of  the  epicyclic  phase  value.  When  the  modal  arm 
amplitudes  are  equal,  these  encounters  equally  retard  and  advance  the  fast  arm  phase  and 
slow  arm  phase,  respectively.  These  card-induced  changes  in  modal  arm  phases  result,  for 
stable  projectile  flight,  in  a  measured  overturning  moment  that  is  larger  than  observed  with 
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spark  ranges. 

For  uniform  yaw  card  arrays,  the  Dirac  delta  approach  is  used  to  develop  difference 
equations  for  a  three-card  sequence.  The  difference  equations  have  as  parameters  the  phase 
interval  between  cards  and  a  value  that  depends  upon  both  the  phase  interval  value  and  the 
stability  parameter.  The  solution  obtained  describes  the  yaw  arm  phases  and  amplitudes 
a.s  the  projectile  passes  through  the  card  array.  It  depends  only  on  the  phase  encounter 
value  for  the  first  card,  the  phase  interval  between  cards,  and  the  stability  parameter  value. 
Stability  boundaries  are  delineated  and  compared  with  approximate  stability  boundaries. 
Conditions  for  maximum  yaw  growth  and  decay  with  each  card  encounter  are  obtained. 

If  the  intervals  between  yaw  cards  are  such  that  the  projectile  experiences  unstable  flight, 
the  encounter  value  for  the  projectile’s  epicyclic  phase  will  asymptotically  approach  a  value 
that  depends  only  on  the  phase  interval  value  between  cards.  This  value  is  called  the  sta¬ 
tionary  encountered  phase  value  for  the  yaw  card  range  and  can  be  readily  calculated.  A 
simple  quaai-universal  relationship  between  the  stationary  phase  value  and  the  card  interval 
is  further  developed.  A  simple  expression  is  also  obtained  for  the  epicyclic  phase  value  that 
occurs  at  each  card  encounter  as  the  projectile  traverses  the  range  for  an  arbitrary  initial 
epicyclic  phase  encounter  value. 
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